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Recent years have seen an explosive research and development of nanoplasmonics in the visible and 
near-infrared (near-ir) frequency regions.— One of the most fundamental effects in nanoplasmonics is 
nano-concentration of optical energy. Plasmonic nanofocusing has been predicted^ and experimen- 
tally achieved . 3 ' 4 ' 5 Nanoconcentration of optical energy at nanoplasmonic probes made possible op- 
tical ultramicroscopy with nanometer-scale resolutions^ and ultrasensitive Raman spectroscopy. 9 
It will be very beneficial for the fundamental science, engineering, environmental, and defense appli- 
cations to be able to nano-concentrate terahertz radiation (frequency 1 — 10 THz or vacuum wave- 
length Ao = 300 — 30 fim). This will allow for the nanoscale spatial resolution for THz imaging 10 
and introduce the THz spectroscopy on the nanoscale, taking full advantage of the rich THz spec- 
tra and submicron to nanoscale structures of many engineering, physical, and biological objects 
of wide interest: electronic components (integrated circuits, etc.), bacteria, their spores, viruses, 
macromolecules, carbon clusters and nanotubes, etc. In this Letter we establish the principal limits 
for the nanoconcentration of the THz radiation in metal/dielectric waveguides and determine their 
optimum shapes required for this nanoconcentration We predict that the adiabatic compression of 
THz radiation from the initial spot size of Rq ~ Ao to the final size of R = 100 — 250 nm can be 
achieved with the THz radiation intensity increased by a factor of xlO to x250. This THz energy 
nanoconcentration will not only improve the spatial resolution and increase the signal/noise ratio 
for the THz imaging and spectroscopy, but in combination with the recently developed sources of 
powerful THz pulses^ will allow the observation of nonlinear THz effects and a carrying out a 
variety of nonlinear spectroscopies (such as two-dimensional spectroscopy), which are highly infor- 
mative. This will find a wide spectrum of applications in science, engineering, biomedical research, 
environmental monitoring, and defense. 



There are existing approaches to deep subwave- 
length THz imaging and probing based on sharp tips 
irradiated by a THz source^ adiabatically-tapered 
metal-dielectric waveguides^ 3 - similar to optical adia- 
batic concentrators*^!^ and nonlinear microscopic THz 
sources^ For the development of the THz nanotechnol- 
ogy, it is extremely important to understand spatial lim- 
its to which the THz radiation energy can be concen- 
trated (nanofocused). 

A major challenge for the nanoconcentration of the 
electromagnetic energy in the THz region is the large ra- 
diation wavelength in vacuum or conventional dielectrics, 
Ao = 30 /im — 300 /im, where the THz radiation can 
only be focused to the relatively very large regions of 
size ~ Ao/2. The developed field of optical energy con- 
centration, which is based on surface plasmon polaritons 
(SPPs), suggests that one of the ways to solving this 
problem is to employ the surface electromagnetic waves 
(SEWs). In the far infrared (ir), the dielectric permit- 
tivity of metals has large imaginary part which dom- 
inates over its negative real part^ This implies that 
SEWs propagating along a metal-dielectric flat interface 
in this frequency range, known as Sommerfeld-Zenneck 
waves ; 16 ' 17 are weakly bound to the surface 18 and can 
hardly be used for the confinement of THz radiation. 

It has been suggested that periodically perforating flat 
surfaces of ideal metals with grooves or holes leads to the 
appearance of SEWs, which mimic ("spoof") SPPs to be 



stronger bound to the surface a 19 ' 20 ' 21 permitting a better 
control over the THz fields. It has been predicted that 
SPPs on an array of parallel grooves cut on the surface 
of a perfect conductor wire can be localized by adiabatic 
deepening of the grooves^ At a point, where grooves are 
approximately a quarter of wavelength of light, the high- 
est concentration is achieved to be on the order of tens 
of micrometers. This method restricts the localization 
point to a particular frequency, making the concentration 
very narrow-band. Also, the depth of a groove should be 
~ Aq/4, i.e., in the tens to hundred micron range, which 
precludes completely nanoscale devices. 

It is well known from microwave technology that the 
ideal-metal waveguides with smooth surfaces support 
TEM waves, where the electric field lines are either in- 
finitely extended or terminate at the metal surfaces nor- 
mally to them. The latter case requires the waveguide 
cross-section topology to be more than single-connected; 
an example may be a coaxial waveguide ("coax"). Such 
waveguides possess are very wide-band in frequency. The 
THz waveguides can be adiabatically tapered to concen- 
trate energy. The idea of adiabatic energy concentration 
comes from ultramicroscop y 23 ' 24 ' 25 and nanoplasmonics, 2 
where it has been developed both theoretically and 
experimentally?^' 5 - and used in ultrasensitive surface en- 
hanced Raman spectroscopy^ Employing these ideas of 
the adiabatic concentration and using a tapered metal- 
dielectric waveguide, the THz spatial resolution achieved 
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is ~ 20 /J,m across the entire THz spectrum*^ 

In this Letter, for the first time, we establish the 
fundamental limits and find the principles of designing 
the optimum and efficient metal/dielectric waveguides 
suitable for the THz nanofocusing. The specific exam- 
ples are for the wide-band concentrators: a plasmonic 
metal wedge cavity and tapered coax waveguides, which 
are terminated by funnel-type adiabatic tapers. Such 
nano-concentrators along with the advent of high-power 
source o 11 ' 26 and sensitive detectors 2 ^ of THz radiation, 
will open up an extremely wide range of possible THz 
applications, in particular, in material diagnostics, probe 
nanoimaging, biomedical applications etc. - cf. Refs. 

Note that an alternative approach to the THz energy 
concentration using doped semiconductor tapers has also 
been proposed.— However, the required heavy doping of 
the semiconductors may cause fast electron relaxation 
due to the collisions with the inflicted lattice defects and 
bring about high losses. Therefore, in this Letter we 
will pursue the adiabatic nanoconcentration of the THz 
radiation using metal/dielectric structures. 

Conventionally for THz and microwave regions, the 
metals are considered as ideal which is equivalent to ne- 
glecting their skin depth l s = Xq /Re y/—e m , where e m is 
the permittivity of the metal (we take into accoint that 
in the THz region |e m | ^> I), and Ao = c/ui is the reduced 
wavelength in vacuum. It is true that in the THz region 
l s = 30 — 60 nm. i.e., l s <C Xq- However, as we show 
below in this Letter, it is the finite skin depth, though 
as small as it is, that principally limits the ultimum lo- 
calization size of the THz fields. For larger waveguides, 
the THz wave energy is localized mostly in the vacuum 
(dielectric) and its losses, which occur in the metal's skin 
layer, are correspondingly small. The effective quality 
factor (or, figure of merit) of the waveguide, which shows 
how many periods the wave can propagate without sig- 
nificantly loosing its energy, can be estimated as 
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where a is the characteristic minimum size of the waveg- 
uide; this estimate becomes a good approximation for a 
metal-dielectric-metal planar waveguide [see below Eq. 
(J2J], When the waveguide size reduces to become on the 
order of the skin depth, a < l s , the THz field is pushed 
into the metal, and the quality factor reduces to Q < 1, 
which implies strong losses. Qualitatively, this estab- 
lishes the limit to the nanoconcentration: for upper THz 
region a > l s « 30 nm, while for the I THz frequency 
a ^ Is ~ 60 nm. These are the practical limits of the THz 
nanoconcentration for the noble metals (silver, gold, and 
platinum) and for aluminum. 

If one pursues the goal of creating enhanced local fields 
in a small region, but not necessarily to efficiently trans- 
fer the THz energy from the far field to the near field, 
then the apertureless SNOM approach, where a sharp 
metal or dielectric tip is irradiated by THz radiation, can, 
in principle, achieve even higher resolution^ However, 
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FIG. 1: Geometry and properties of the THz TM mode in 
a parallel metal-slab waveguide, (a) Schematic of the waveg- 
uide. The width of the dielectric gap o and the skin depth l s 
are indicated, (b) An instantaneous distribution of the lon- 
gitudinal electric field E y along the propagation coordinate y 
for a — 10 fim and frequency 1 THz in a silver-vacuum-silver 
waveguide, (c) The same as in panel (b) but for a = 200 
nm. (d) Modal refraction index n — k/ko (Ren is denoted by 
the red line and Imn by the blue line) as a function of the 
waveguide width a. Dashed green line indicates the value of 
n for the perfect conductor. Skin-depth value is shown by the 
vertical dashed line. 



the efficiency of utilizing the THz energy of the source in 
this case will be extremely low; the stray, far-field THz 
energy may create a significant parasitic background. 

Here and below in this Letter, we consider examples 
of the THz adiabatic nanoconcentration quantitatively, 
where the effect of the specific geometry will become ap- 
parent. Consider first a parallel plate waveguide that 
consists of a dielectric slab of thickness a with dielec- 
tric permittivity Ed sandwiched between two thick metal 
plates (with thickness of at least a few l s , i.e., greater 
than 200 nm in practical terms) [see FigQJa)]. The per- 
mittivity of the metal e m in the THz region has a very 
large (> I0 6 ) imaginary part that defines the very small 
skin depth l s < 100 nm, which justifies the usual consid- 
eration of the metals as perfect conductors^ However, 
as we have already mentioned, for our purposes of the 
THz nanoconcentration, we need to take into account 
the field penetration into the metal, which makes the 
problem plasmonic. In this case, the propagating modes 
of the system are SPPs, which are TM modes character- 
ized by the symmetry with respect to the reflection in 
the center plane. We will orient the coordinate system 
with its z axis normal to the plane and the y axis in the 
direction of propagation. The symmetric (even) modes 
have even field components H x and E z and odd E y ; the 
parity of the antisymmetric (odd) modes is opposite. 

From plasmonics it is known that the even modes have 
a larger fraction of their energy localized in the dielec- 
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trie and the odd modes in the metal. Therefore, the 
even modes have much smaller damping and are, there- 
fore, most suitable for the THz energy concentrations. 
The dispersion relation for the even modes is given in 
the Methods section as Eq. ([5]). This equation can be 
much simplified and solved in a closed analytical form 
taking into account that we are interested in the sub- 
wavelength focusing, i.e., a < A, where A = Xo/ed is the 
reduced wavelength in the dielectric; also, in the entire 
THz region l s <C A. This shows that there exists a small 
parameter in the problem l s a/X 2 <C 1 [see also Eq. @], 
which allows one to solve analytically the dispersion rela- 
tion (0 obtaining the modal refraction index n = k/ko, 
where k is the THz wavevector, and k a = 1/Aq, 



the Methods section. The WKB approximation is appli- 
cable under the conditions that 
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where the approximate equality is valid for not too tight 
nanofocusing, i.e., for l s <C a. From this, we can obtain 
the quality factor of the waveguide Q = Ren/Imn = 
2a/ l s , giving a quantitative meaning to the estimate (fTj). 

Plasmonic effects (i.e. those of the finite skin depth) 
are illustrated in Fig. [T] for silver-vacuum-silver waveg- 
uide and frequency of 1 THz. Panels (b) and (c) display 
the longitudinal electric field E y obtained by the exact 
solution of the Maxwell equations. Note that this field 
component is absent for the ideal conductor; here it is 
relatively small: on the order of 10~ 3 of the transverse 
field. Panel (b) illustrates the case of a relatively wide 
waveguide (a = 10 /an), where it is evident that the elec- 
tric field is localized mostly in the dielectric region of the 
waveguide, and the extinction of the wave is small. In a 
sharp contrast, for a nanoscopic waveguide (a = 200 run) 
in panel (c), the electric field significantly penetrates the 
metal. In accord with our arguments, there is a very 
significant extinction of the fields as they propagate; the 
retardation effects are also evident: the lines of equal 
amplitude are at an angle relative to the normal (z) di- 
rection. The dependence of the modal refraction index 
on the thickness a of the waveguide obtained from Eq. 
@ is plotted in Fig.[T](d). This index increases as a be- 
comes comparable with the skin depth. While Ren and 
Imn increase by the same absolute amount, the quality 
factor Q, obviously, greatly decreases with decrease of a. 
The mode described above can be used for broadband 
energy concentration of THz waves. 

To introduce the THz nanoconcentration, consider a 
metal-dielectric-metal waveguide that is slowly (adiabat- 
ically) tapered off as a wedge, as illustrated in Fig. [3] 
(a). Because of the adiabatic change of the parameters, 
a wave propagating in such a waveguide will adjust to 
it without reflection or scattering, just as it takes place 
in nanoplasmonic waveguides^ As a result, propagating 
it will concentrate its energy, conforming to the taper- 
ing of the waveguide. The corresponding solution can be 
obtained from the Maxwell equations using the Wentzel- 
Kramers-Brillouin (WKB) approximation, similarly to 
the nanoplasmonic case in the visible^ as described in 



6=\d (Refc -1 ) /dy\ < 1 , \da/dy\ < 1 
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where S is the well-known adiabatic parameter describing 
how slowly the modal wavelength changes on a distance 
of its own, and \da/dy\ is a parameter describing how 
adiabatically the transverse size of the confined mode 
changes along the propagation coordinate. 

In the WKB approximation, the behavior of the domi- 
nating transverse field component E z as a function of the 
coordinate y along the propagation direction is shown for 
the last 6 /jm of the propagation toward the edge in Fig. 
[5](b). There is a clearly seen spatial concentration of the 
energy and increase of the field as the wave is guided into 
the taper. The predicted behavior of the two components 
of electric field and the magnetic field for the last micron 
of the propagation is shown in panels (c)-(e). It appar- 
ently indicates the adiabatic concentration, without an 
appreciable loss of the intensity. The THz wave follows 
the waveguide up to the nanometric size. 

The red line in Fig. [3] (f ) indicates that the local inten- 
sity I as the function of the thickness a of the waveguide 
for a < 4 /im increases significantly with l/a, in qual- 
itative accord with the behavior expected for the neg- 
ligibly low losses. This intensity reaches its maximum 
for a — 1.6 fim and then starts to decrease as the losses 
overcome the adiabatic concentration. At smaller thick- 
nesses, a < 400 nm, the intensity in Fig. [2] (f) starts 
to increase again, which is unphysical. The reason is 
revealed by the behavior of the adiabatic parameter S 
shown by the blue line: for a < 400 nm, S becomes rel- 
atively large (comparable with 1), i.e., the adiabaticity 
is violated. This is due to the fact that the fraction of 
the THz field energy propagating in the metal is dramat- 
ically increased for a < 400 nm due to the constricted 
transverse extension of the dielectric in the waveguide. 
This causes a significant loss per wavelength A, leading 
to a rapid change of the wave vector fc, breaking down 
the adiabaticity. This constitutes a fundamental differ- 
ence from the nanoplasmonic adiabatic concentration in 
the optical region where the adiabatic parameter is con- 
stant, and the adiabaticity holds everywhere including 
the vicinity of the tip^ 

To provide for the optimum guiding of the THz wave 
and its concentration on the nanoscale, the terminat- 
ing (nanoscopic) part of the waveguide should be ta- 
pered slower, in a funnel-like manner. That is, one 
needs to decrease the grading da/dy of the waveguide 
near the edge in order to keep the adiabaticity parame- 
ter 5 = Id (Refc -1 ) I da x da/dy\ approximately constant 
and small enough to prevent the back-reflection. Be- 
cause for the adiabatic grading (tapering) , the derivative 
d (Refc -1 ) I da does not depend on the grading (it is the 
same as for the plane waveguide) and is only a function 
of a, the equation 5 = S(y) is a differential equation for 
the shape of the waveguide that can be easily integrated. 
This results in the dependence of the thickness a on the 
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FIG. 2: Adiabatic concentration of THz field energy in a 
graded waveguide, where a dielectric wedge is surrounded by 
the thick silver layer, (a) Schematic of energy concentration, 
where 9 is the wedge opening angle, the arrow indicates the 
direction of propagation of the THz wave, and the red high- 
lights the area of the adiabatic concentration. The orienta- 
tion of the coordinate system is shown in the inset, (b) An 
instantaneous distribution of the transverse electric field E z 
of the THz wave propagating and concentrating along the 
wedge waveguide for the last 6 mm of the propagation to- 
ward the edge. Note the difference in scales in the z and y 
directions, (c) An instantaneous spatial distribution of the 
transverse electric field E z close to the edge of the wedge, 
for the last 640 /im of the propagation, (d) The same as (c) 
but for the longitudinal (with respect to the propagation di- 
rection) component of the field E y . (e) The same as (c) but 
for the transverse component of the magnetic field H x . The 
units of these field components are arbitrary but consistent 
between the panels, (f) Dependence of THz field intensity in 
the middle of waveguide on the dielectric gap width a (the 
red line) . The blue curve displays the dependence on a of the 
adiabatic parameter S, scaled by a factor of 5. The values of 
a indicated at the successive horizontal axis ticks differ by a 
factor of 10" 1/2 , i.e., by 5 dB. 



longitudinal coordinate y determined by a simple integral 
Ren- 1 (a)=k [ S(y)dy , (4) 



where n(a) is the modal index denned in this case by Eq. 
([2]), and S(y) is the desired dependence of the adiabatic 
parameter along the waveguide, which is an arbitrary 
function of y satisfying the adiabaticity conditions ©. 
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FIG. 3: Terahertz energy concentration in adiabatically ta- 
pered curved-wedge waveguide, (a) Instantaneous distribu- 
tion of the transverse component of the THz electric field E z 
(in the central plane z — 0) as a function of the coordinate 
y along the propagation direction for the last 400 /im of the 
propagation, (b) The same as in panel (a) but for the longi- 
tudinal electric field component E y . (c) The same as panel 
(a) but for the transverse magnetic field H x . The units for 
the fields are arbitrary but consistent between the panels, (d) 
The THz field intensity I (relative to the intensity 7o at the 
entrance of the waveguide) as a function of the dielectric gap 
thickness a is shown by the red line. The adiabatic parame- 
ter scaled by a factor of 10 as a function of a is indicated by 
the blue line. The values of a indicated at the horizontal axis 
ticks correspond to the values of y at the ticks of panels (a)- 



The geometry of an adiabatically-tapered end of the 
silver/vacuum waveguide found from Eq. ([¥]) and satis- 
fying Eq. ([3]) and the corresponding WKB solutions for 
the 1 THz fields are shown in Fig. [3] (a)-(c). The opti- 
mum shape of the waveguide in this case is funnel-like, 
greatly elongated toward the edge. The nanoconcentra- 
tion of the field is evident on panels (a)-(c), as well as 
its penetration into the metal for a < 100 nm. As these 
panels show quantitatively and the red curve on panel 
(d) qualitatively, the field intensity reaches its maximum 
at a sa 300 nm where it is enhanced with respect to the 
field at the entrance to the funnel waveguide by a modest 
factor of 1.2. At the same time, the adiabatic parameter 
5 decreases toward the tip from 0.07 to 0.05, indicated 
the applicability of the WKB approximation everywhere. 
Note that the this funnel-shaped wedge, indeed, contin- 
ues the linearly-graded wedge waveguide shown in Fig. [21 
which yields the enhancement factor of ss 8 at a = 2 /an. 
Sequentially, these two waveguides provide the intensity 
enhancement by approximately x 10 while compressing 
the THz wave to the thickness of a = 300 nm and the 
enhancement by a factor of 3 for a = 100 nm. 

Thus, true nanolocalization of THz radiation in one 
dimension (Id) is possible. The minimum transverse size 
of this nanolocalization is determined by the skin depth, 
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as we have already discussed qualitatively in the intro- 
ductory part of this Letter. The obtained Id beam of the 
nanoconcentrated THz radiation may be used for differ- 
ent purposes, in particular as a source for the diffraction 
elements including the nanofocusing zone plates of the 
type introduced in Ref. I33L 

The two-dimensional (2d) concentration of the THz ra- 
diation can be achieved by using an adiabatically-tapered 
conical coax waveguide, whose geometry is illustrated in 
Fig. 0] (a). The central metal wire of radius r is sur- 
rounded by a dielectric gap of the radial thickness a, 
which is enclosed by a thick (~ 200 nm or thicker) outer 
metal shell. Both r and a are smooth functions of the 
longitudinal coordinate y, which describes the tapering 
of the coax toward the apex (tip) at y = 0. The THz 
waves propagate from the wide end of the coax toward 
the apex, adiabatically following the tapering. In the 
spirit of WKB, for any particular y the wave behavior 
for the tapered coax is the same as for a cylindrical coax 
with the values of r and a equal to the local values r(y) 
and a(y). 

The dispersion relation for the coax waveguide that 
takes into account the plasmonic effects (i.e., the pen- 
etration of radiation into the metal and the concurrent 
losses) is obtained in the Methods section as Eq. ([9]). 
Calculated from this expression, the dependence of the 
modal refractive index n = k(a,r)/ko on the dielectric 
gap a is displayed in Fig. 0] (b) for the frequency of 1 
THz, silver as a metal, and vacuum in the dielectric gap. 
The results are shown for two values of the radius of the 
central wire: r — 10 /jm and r = 60 nm. As one can 
see, the real part of the modal index practically does not 
depend on r; it starts growing when a decreases. The 
imaginary part of the index n increases when the central 
wire thickness r decreases, but this dependence is very 
weak. Both Ren and Imn grow dramatically for r < l s . 
This is due to the penetration of the THz field into the 
metal, i.e., it is a plasmonic effect. 

The WKB solution for the radial field E p in the cross 
section of this coax waveguide is shown for the last 3 jirsx 
of the propagation toward the tip in Fig. Q] (c). The adi- 
abatic following and energy concentration are evident in 
this panel. The penetration into the metal of the tan- 
gential (to the metal surface) field components E y and 
H v is noticeable in Figs. 2] (d) and (e). The intensity / 
of the THz field (relative to the intensity Iq at the en- 
trance of the waveguide) as a function of the waveguide 
outer radius R = r + a is shown by the red line in Fig. 0] 
(f). Dramatically, it shows the adiabatic nanoconcentra- 
tion and the intensity increase by more than two orders 
of magnitude for the nanoconcentration from the waveg- 
uide radius R — 300 /im, where the THz radiation can 
be focused, toward R — 300 nm. However, the dramatic 
increase of the adiabatic parameter 5 [plotted by the blue 
line in Fig. [4] (f )] for R < 1 /Ltm shows that these results 
can only be trusted for R > 1 /im. 

For the true 2d nanoconcentration of the THz radia- 
tion below this micron-scale radius, similar to the Id case 
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FIG. 4: Geometry, modal index of refraction, and THz energy 
concentration in conically-tapered metal-dielectric waveguide, 
(a) Schematic of geometry and energy concentration. The 
central wire and the coax shell are shown along with the 
schematic of the THz energy concentration, (b) Dependence 
of modal refraction index n in coaxial waveguide on the di- 
electric gap width a for two central wire radii: r = 10 /im 
and r — 60 nm. The color coding of the lines is indicated. 
The dielectric in the gap is vacuum, (c) Instantaneous dis- 
tribution of the radial THz electric field amplitude E p in the 
cross section of the coax for the last 3 mm of the propagation 
toward the tip. The amplitude of the field is color coded by 
the bar at the top of the panel, (d) Instantaneous distribu- 
tion of the longitudinal THz electric field amplitude E y on the 
coordinate y for the last 620 /im of the propagation, (e) The 
same as (d) but for the transverse magnetic field H v . The 
units of these field components are arbitrary but consistent 
between the panels, (f) Dependence of THz field intensity in 
the middle of waveguide gap on the waveguide outer radius 
R = r + a is shown in red. The blue curve displays the adi- 
abatic parameter S as a function of R, scaled by a factor of 
10 2 . The values of R indicated at the successive horizontal 
axis ticks differ by a factor of 10 _1//2 , i.e., by 5 dB. 



of the wedge, to preserve the adiabaticity, a funnel-like 
tapering is necessary. Generally, the tapering of the cen- 
tral wire and that of the outer metal shell do not need 
to be the same. However, we found that better results 
are obtained when it is the case, i.e., the waveguide is 
tapered-off self-similarly. In specific calculations, as ev- 
erywhere in this Letter, we assume that the metal of 
the waveguide is silver, the dielectric is vacuum, and the 
frequency is 1 THz. Doing so, we have found the cor- 



6 




600 400 200 1000 400 220 150)10 85 65 



r(^im) /f(nm) 

FIG. 5: Adiabatic terahertz energy concentration in a self- 
similarly curved, funnel-shaped coaxial waveguide, where the 
metal is silver, and the dielectric in the gap is vacuum. The 
dielectric gap is between the pair of the neighboring curved 
lines, and the metal is everywhere else, (a) Instantaneous 
distribution of the radial (transverse) component E p of the 
electric field of the guided THz wave as a function of the 
propagation coordinate along the wedge y for the last 600 /im 
of the propagation, (b) The same for the longitudinal electric 
field component E v . (c) The same for transverse magnetic 
field H Vl whose lines form circles around the central metal 
wire. The units of these field components are arbitrary but 
consistent between the panels, (d) The THz intensity J as a 
function of the waveguide radius R, displayed relative to the 
intensity 7o at the beginning of the waveguide (red line) . Adi- 
abatic parameter S multiplied by a factor of 10 as a function 
of R (blue line) . The values of the radius R shown at the ticks 
correspond to those of y shown in panels (a)-(c). 

responding grading of the waveguide using Eq. (J3|) and 
setting S — 0.05, which is small enough to satisfy the 
adiabaticity very well. In this case, indeed, we have used 
the corresponding dispersion relation The obtained 
shape of the waveguide is a strongly-elongated funnel, as 
shown in Fig. [5] (a)-(c). These figures display the THz 
fields that we have calculated in the WKB approximation 
for this waveguide. As one can see from these figures, 
within the last half micron of the propagation, the elec- 
tric and magnetic fields of the THz wave efficiently fol- 
low the adiabatically curved waveguide. The penetration 
into the metal of the tangential (to the metal-dielectric 
interfaces) field components for y < 400 /im is evident in 
panels (b) and (c). The longitudinal electric field com- 
ponent E y is significantly localized in the central metal 
wire [panel(b)], which is a plasmonic effect. 

The dependence of the THz field intensity in the gap 
(relative to the intensity Iq at the entrance of this fun- 
nel) on the total radius of the waveguide R is shown in 
Fig. [5] (d) by the red line. In this case, the adiabatic 
concentration is very efficient. The intensity of the THz 
radiation increases by a factor of x5 when it is com- 
pressed from the initial radius of R = 1 /im to the radius 



R » 250 nm. The penetration of the fields into the metal 
for smaller values of the radius R (tighter confinement) 
causes losses that dominate over the effect of the concen- 
tration. Again, we remind that this funnel waveguide is 
a continuation and termination for the straigt cone that 
yields the field enhancement by x50 for R = 1 /im [see 
Fig. Q] (f )] . Consecutively, these two waveguides (the ini- 
tial cone continued and terminated by the funnel) are 
very efficient, adiabatically compressing the THz radia- 
tion from the initial radius R = 300 fim to the radius 
R = 250 nm increasing its intensity by a factor x250. 
Even for the final radius R = 100 nm, the total THz 
intensity is increased by a factor of xlO (which is the 
products of factors x50 for the cone part and x0.2 for 
the funnel. Thus, the optimally graded plasmonic-metal 
2d waveguide is very efficient in the concentration and 
guidance of the THz fields with the transverse radius of 
confinement R < 100 nm. 

To discuss the results, we have shown that the THz ra- 
diation can be concentrated to the ~ 100 nm transverse 
size in adiabatically graded plasmonic (metal/dielectric) 
waveguides. In the optimum adiabatically-graded, coax- 
ial waveguide, which consists of the initial cone termi- 
nated with a funnel, the radiation of a 1 THz frequency 
whose wavelength is 300 /im, can be compressed to a 
spot of 250 nm radius, where its intensity increases by a 
factor of x250. Even in the case of the extreme compres- 
sion to a spot of the 100 nm radius, the THz intensity is 
enhanced by one order of magnitude with respect to the 
initial intensity of the 300 /im spot at the entrance of the 
waveguide. The physical process that limits the extent of 
this spatial concentration is the skin effect, i.e., penetra- 
tion of the radiation into the metal that causes the losses: 
the THz field penetrates the depth of l s = 30 — 60 nm of 
the metal, which by the order of magnitude determines 
the ultimum localization radius. 

The THz nanoconcentration predicted in this Let- 
ter for optimally-graded adiabatic plasmonic waveguides 
provides unique opportunities for THz science and tech- 
nology, of which we will mention below just a few. The 
nanoconcentration of the THz radiation will provide the 
THz ultramicroscopy with a THz source of unprece- 
dented spatial resolution and brightness. The increase 
of the THz intensity by two orders of magnitude along 
with the novel high-power THz sources^ would allow 
the observation of a wide range of electronic and vibra- 
tional nonlinear effects in metal, semiconductors, insula- 
tors, and molecules. 

These nonlinear THz phenomena can be used to inves- 
tigate behavior of various materials in ultrastrong fields, 
for nonlinear spectroscopy (including the multidimen- 
sional spectroscopy), and for monitoring and detection of 
various environmental, biological, and chemical objects 
and threats such as single bacterial spores and viruses. 
Such applications will certainly be helped by very large 
absorption cross sections of various materials in the THz 
region. A distinct and significant advantage of the adia- 
batic nanofocusing is that the THz energy is mostly con- 
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centrated in the hollow region of the waveguide, whose 
size can be made comparable with the size of the objects 
of interest: in the range from 1 micron to 70 nm, which is 
a typical range for bacteria and their spores, and viruses. 
This will assure high sensitivity and low background for 
the objects that are confined inside these waveguides. 

Consider as a specific example the spectroscopy or de- 
tection of single particles, such as, e.g., anthrax spores, in 
the air. A sample containing the suspected nanoparticles 
in a gas, which can be air for the frequencies in the trans- 
parency windows, can be pumped through a THz waveg- 
uide, and the detection can be made for each particle in 
the gas separately on the basis of the two-dimensional 
nonlinear THz spectra that are expected to be highly in- 
formative for the detection and elimination of the false- 
positive alarms. Likewise, many other scientific, tech- 
nological, environmental, and defense applications may 
become possible. 



I. METHODS 

Terahertz TM Wave in Finite-Conductivity 
Parallel Plate Waveguide 



the dielectric gap width. The characteristic relation for 
the TM modes of this waveguide has has the following 
form 



/ Ig(n d 
( K (k, 



r) 



C h(K d r) \ ( 



K ( Km R) 



K x (K d R) \ _ 



K 1 (K m R) 

, p Kx{K, d r) \ I I (n d R) , c I x (K d R) \ 
+ ? /i(K m r) J \K ( Km R.) ^^K^KmR) ) 



(7) 



where I v {x) and K v (x) are modified Bcsscl functions, 



and £ 



This equation is quadratic with respect 



to £ and can be written in the form a£ 2 + /3£ + 7 = 0, 
where the coefficients a, /3, and 7 can be easily found 
by comparison to Eq. ([7|) as combinations of the Bessel 
functions. It can obviously be resolved for £ yielding 



(3 ± V/? 2 - 4a 7 



2a 



(8) 



In the THz region, only the mode with the minus sign 
in Eq. ((5]) propagates. It can be treated in a manner 
similar to the mode described by Eq. (f5]). The equation 
(|8|) can be expanded over the small parameter n^a <C 1, 
and the explicit form of the modal refraction index can 
be readily obtained as 



A parallel plate waveguide supports an even TM mode 
with wavenumber k, which satisfies the dispersion rela- 
tion 

> 2 / £m K d 



tanh i 



(5) 



where Kd = (k 2 — 
terahertz range, e 
1. Therefore, n m 
l/ReK m = %/2/(fco\/k 



edkl) 



— [k Smk^ j 



In the 

is mainly imaginary, where Im e m S> 
koV- £ m = Ij 1 ^ ~ i), where / s = 
is the metal skin depth, which 
is on the order of tens of nanometers. We also assume 
that Kda <C 1, which is always the case for the mode 
under consideration because either this mode is close to 
the TEM mode where k = kgSd, or the gap a is thin 
enough. This leads to a closed expression for the index 

of refraction of the mode, which is Eq. ([2]). Using this, 

1/2 

one can check that k^o, \Edal s /% ) . Consequently, 
the applicability condition of the approximation used is 



(e d aZ s /A 2 ) 1/2 «l 



(6) 



This condition is satisfied for the realistic parameters of 
the problem. For instance, for the frequency / = 1 THz, 
the skin depth for metals is l s ~ 60 nm, while reduced 
wavelength is X = 75 /im. The condition ((6]) is well sat- 
isfied for a <C 100 /im, i.e., in the entire range of interest 
to us. 



Terahertz TM Wave in Finite-Conductivity 
Coaxial Waveguide 



n = Jel 1 



/o(Kmr) , K (n m R) 
/iM + K x {K m R) 



2a 



1/2 



(9) 



where K m = l s 1 {1 — i). Similar to the wedge waveguide 
case, the applicability condition of this solution is given 
by Eq. ©. 



C. WKB Solution for the TM Wave in Graded 
Waveguide 

In the WKB approximation, a solution of the Maxwell 
equations can be represented as a wave with amplitude 
and phase that are slowly varying functions of y on the 
scale of local wavelength. The behavior in the transverse 
direction z is the same as for the non-graded system. The 
WKB solution is valid if the adiabatic parameter is small 



1 



dy k(y) 



< 1 . 



(10) 



In the WKB approximation, the phase of the mode 
(eikonal) is given by an integral 



<K?/) = k n(y)dy 



(11) 



where n(y) is the local refraction index of the mode. The 
behavior of wave amplitude as a function of the propa- 
gation coordinate y is found from the condition of flux 
conservation: 



Consider a coaxial waveguide (coax) with the inner 
wire radius r and the outer radius R = r + a, where a is 



v g(y) / W(y, z)dz = const 



(12) 
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where v g (y) = duj/dk is the wave local group velocity, 
and W(y, z) is energy density in the mode. 
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